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STELLAR VELOCITY DISPERSION AND BLACK HOLE MASS IN THE BLAZAR MARKARIAN 501 

Aaron J. Barth'-^, Luis C. Ho^ and Wallace L. W. Sargent' 



O 
O 
(N 



> 
O 

o 

O 



Oh- 
I 

O 



X 



ABSTRACT 

The recently discovered correlation between black hole mass and stellar velocity dispersion provides a new 
method to determine the masses of black holes in active galaxies. We have obtained optical spectra of Markarian 
501, a nearby 7-ray blazar with emission extending to TeV energies. The stellar velocity dispersion of the host 
galaxy, measured from the calcium triplet lines in a 2" x 3"7 aperture, is 372 ± 18 km s"'. If Mrk 501 follows 
the M, - cr* correlation defined for local galaxies, then its central black hole has a mass of (0.9-3.4) x 10^ Mq. 
This is significantly larger than some previous estimates for the central mass in Mrk 501 that have been based on 
models for its nonthermal emission. The host galaxy luminosity implies a black hole of 6 x 10^ Mq, but this is 
not in severe conflict with the mass derived from a^, because the M, -Lbuige correlation has a large intrinsic scatter. 
Using the emission-line luminosity to estimate the bolometric luminosity of the central engine, we find that Mrk 
501 radiates at an extremely sub-Eddington level of L/LEdd ~ 10"^. Further applications of the M, - cr* relation to 
radio-loud active galactic nuclei may be useful for interpreting unified models and understanding the relationship 
between radio galaxies and BL Lac objects. 

Subject headings: BL Lacertae objects: individual (Mrk 501) — galaxies: active — galaxies: elliptical and 
lenticular — galaxies: kinematics and dynamics — galaxies: nuclei 



1. INTRODUCTION 

The tight correlation recently discovered between stellar ve- 
locity dispersion and black hole mass in nearby galaxy bulges 
(Ferrarese & Merritt 2000; Gebhardt et al. 2000a) has become 
the key to our understanding of black hole demographics, as 
well as a new tool for probing the evolution of galaxies and 
quasars. An equally important aspect of the M, - a^, correla- 
tion is its predictive power While dynamical mass measure- 
ments are observationally difficult and only possible for a lim- 
ited number of galaxies, the M, - cr* relation makes it possible 
to obtain a black hole mass estimate good to ^ 40% accuracy or 
better from a single measurement of bulge velocity dispersion. 

Having such a straightforward method to estimate black hole 
masses is a tremendous boon to studies of active galactic nu- 
clei (AGNs), because M, is a fundamental parameter affecting 
the energetics and emission properties of AGNs. Gebhardt et 
al. (2000b) and Ferrarese et al. (2001) have recently shown that 
black hole masses derived for Seyfert galaxies via the M, - (t^ 
relation are consistent with masses determined by reverbera- 
tion mapping, providing added confidence that AGNs do fol- 
low the same correlation as inactive galaxies. For some classes 
of AGNs, using the correlations between M, and host galaxy 
properties may be the only reliable way to estimate the central 
masses. BL Lac objects fall in this category, since more direct 
methods for determining M, (stellar dynamics or reverberation 
mapping) cannot be applied. 

Markarian 501 is one of the nearest known BL Lac objects. 
Its redshift is z = 0.0337 (Ulrich et al. 1975), corresponding to 
a distance of 144 Mpc for Hq = 70 km s"' Mpc"'. The host 
of Mrk 501, also known as UGC 10599, is a giant elliptical 
that appears morphologically normal in ground-based images 
(Hickson et al. 1982; Abraham, McHardy, & Crawford 1991; 
Stickel, Fried, & Kuhr 1993; Wurtz, Stocke, & Yee 1996; Nils- 
son et al. 1999). Stellar absorption lines are visible in the nu- 
clear spectrum, along with weak emission lines (Ulrich et al. 
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1975; Stickel et al.l993). Mrk 501 is also a 7-ray source, and is 
one of the few extragalactic objects from which TeV emission 
has been detected (Quinn et al. 1996; Bradbury et al. 1997). 

In this Letter, we present new optical spectra of Mrk 501. 
The Ca II triplet lines are clearly detected, albeit substantially 
diluted by the nonthermal continuum. The stellar velocity dis- 
persion measured from these lines is 372 ± 18 km s"'. Using 
this result in conjunction with the M, - a^, correlation, we de- 
rive an estimate of the mass of the black hole in Mrk 501. 



2. OBSERVATIONS AND REDUCTIONS 

Observations were obtained with the Palomar Hale 200-inch 
telescope and Double Spectrograph (Oke & Gunn 1982) on 
2001 June 24 UT. On the red side, we used a 1200 grooves 
mm"' grating blazed at 9400 A, and a wavelength setting giv- 
ing coverage of 8435-9075 A at 0.635 A pixel"'. For the blue 
side, we used a 600 grooves mm"' grating blazed at 4000 A, 
covering 4200-5950 A at 1.72 A pixel"'. A 2"-wide sHt was 
used for all observations. The total exposure time for Mrk 501 
was 3.5 hours, broken into individual 30-minute exposures. 

The airmass ranged between 1.01 and 1.20 during the ob- 
servations, and the slit was held at a fixed position angle of 
49°. There were intermittent thin clouds during the observa- 
tions, and the seeing was typically ^ l."5. We observed Mrk 
501 again on June 26, but due to thick clouds the data were of 
very poor quality and we do not include these exposures in our 
analysis. We also observed 24 stars as velocity templates dur- 
ing the observing run. The stars included giants, subgiants, and 
dwarfs with spectral types between F3 and K5. 

Each frame was bias-subtracted and flat-fielded using stan- 
dard techniques. Spectral extractions were performed using a 
width of 3 "7, corresponding to 8 pixels on the red side and 6 
pixels on the blue side. The extractions were flux-calibrated 
using observations of the standard star BD H- 174708, and 
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wavelength-calibrated using exposures of a He-Ne-Ar lamp, 
with a final linear shift based on the wavelengths of night sky 
emission lines. Telluric absorption bands longward of 8925 
A were removed by dividing by the normaUzed spectrum of 
the standard star, but the corrected spectra are extremely noisy 
longward of 8950 A. A few narrow, weak telluric features with 
equivalent widths of < 0.04 A remain in the object spectra in 
the region 8800-8925 A. These features could not be removed 
by standard star division since the standard star has intrinsic ab- 
sorption features in this wavelength range that are comparable 
to, or stronger than, the telluric features. 

The spectrum of Mrk 501 is displayed in Figure 1, along with 
the nearby elliptical galaxy NGC 4278 as a comparison. Sev- 
eral stellar absorption features are clearly visible in the Mrk 501 
spectrum, despite the strong dilution by nonthermal emission. 
These include the Ca II triplet lines at AA8498, 8542, 8662 A, 
and on the blue side, the G band and Mg lb. Weak emission 
lines of |0 III| AA4959,5007 are also detected in Mrk 501, but 
H/3 emission is absent. The Ca II A8542 absorption line has an 
equivalent width of only ~ 1 A (in comparison with 3-4 A in 
our K-giant stellar templates), indicating that the spectrum is 
dominated by nonthermal emission. 

3. MEASUREMENT OF VELOCITY DISPERSION 

To measure the velocity dispersion, we used a template fitting 
routine based on the method described by van der Marel (1994), 
which performs a direct fit of broadened stellar templates to 
the galaxy spectrum. The object and template spectra are first 
rebinned to a wavelength scale that is linear in log A. The tem- 
plate is broadened by convolution with a Gaussian, and its spec- 
tral shape is adjusted with an additive dilution, representing the 
nonstellar continuum, and multiplication by a quadratic poly- 
nomial, to allow the continuum shape of the template to match 
the object. The galaxy's radial velocity is an additional free pa- 
rameter in the fit. This procedure is repeated using a downhill 
simplex search algorithm (Press et al. 1988), allowing all pa- 
rameters to vary freely, to find the parameter values that mini- 
mize x^, yielding the best-fitting Gaussian width a-^. The fitting 
routine was tested with spectra of velocity dispersion standard 
galaxies from the catalog of McElroy (1995) that were observed 
during the same run, and in all cases our results agreed well 
with the previously measured values. 

The assumption of a Gaussian for the hne-of-sight veloc- 
ity distribution is a reasonable choice given the relatively low 
signal-to-noise ratio of the Mrk 501 spectrum; our data do not 
have the sensitivity needed to detect higher-order moments of 
the velocity distribution. The high-velocity wings on the veloc- 
ity profiles due to the central black hole should not be signifi- 
cant given the large spectroscopic aperture, which corresponds 
to 1.3 kpc X 2.0 kpc at the distance of Mrk 501. 

For the red side spectrum, we chose a fitting region just large 
enough to contain the Ca lines, 8475-8680 A. The spectrum 
becomes extremely noisy beyond 8680 A and we truncated the 
fit at that wavelength; the red wing of the A8662 line is es- 
sentially lost in this portion of the spectrum due to the telluric 
water absorption band, even after division by the standard star. 
We excluded a few remaining weak telluric features from the fit, 
including one located in the blue wing of the A8542 absorption 
feature. The best-fitting template is HD 188056, a K3m star 
(Figure 2). With this template, the derived velocity dispersion 
is cr* = 372 ± 7 km s~\ and the standard deviation among all 
templates observed during the run is 17 km s~^ Since template 



mismatch is the dominant source of uncertainty in the measure- 
ment, we combine the template-matching and fitting uncertain- 
ties in quadrature for a final result of = 372 ±18 km s~' . 

The Ca triplet region is expected to be less sensitive to stel- 
lar population variations than the blue spectral region (e.g.. 
Dressier 1984), so the fit to the Ca lines should give the best 
estimate of the velocity dispersion. We also performed fits to 
several portions of the blue spectrum, and found that the best 
results were obtained for a region just redward of the Mg lb 
line. The Mg lb line itself did not match the template spectra 
well, however, and it was excluded from the final fits. Over the 
range 5200-5600 A, the best-fitting template was the K3III star 
HD 125560, giving cr* = 386 ± 9 km s~\ in reasonable agree- 
ment with the red side results. However, several of the template 
stars gave very poor fits, and the standard deviation among aU 
templates for this region was 81 km s"^ 

4. THE MASS OF THE BLACK HOLE IN MRK 501 

If we assume that Mrk 501 follows the M, - relation de- 
fined for nearby galaxies, then it is straightforward to convert 
our measurement of cr* into an estimate of M, . This assumption 
is supported by the fact that the Mrk 501 host galaxy appears 
morphologically undisturbed, without any indications of recent 
major merger or interaction events that might cause the galaxy 
to deviate from the correlation. 

The most recent fits to the M, - relation are given by Ko- 
rmendy & Gebhart (2001), who find M, cx a^^^, and by Mer- 
ritt & Ferrarese (2001), who find M, oc cr^^^. Given the dis- 
agreement between these two groups, the most conservative 
approach is to apply both relations and derive a range of pos- 
sible black hole masses for Mrk 501. Another issue is the fact 
that Gebhardt et al. (2000a) use the luminosity-weighted ve- 
locity dispersion within the half-light radius, while Ferrarese & 
Merritt (2000) base their fit on the "central" velocity dispersion 
measured within an aperture of radius r^/S. For our purposes, 
the distinction between these two choices is not large; Gebhardt 
et al. (2000a) show that variations in luminosity-weighted ve- 
locity dispersion for apertures of different size rarely exceed 
10% for r < 4re. One possibiUty would be to use the pre- 
scription of J0rgensen, Franx, & Kjaergaard (1995) to normalize 
our measured aperture dispersion to or r<,/8. However, this 
would introduce additional uncertainty as the value of r,, for 
Mrk 501 is not well constrained by existing ground-based imag- 
ing data. Values for re compiled from the Uterature by Nilsson 
et al. (1999) range from 9" to 20". Thus, we use our measured 
velocity dispersion without applying any aperture corrections. 

Propagating the uncertainty in a^, as well as the uncertain- 
ties in the coefficients of the M, - relation, the black hole in 
Mrk 501 has M, = (1 .3 ± 0.4) x 10^ Mg with the Kormendy & 
Gebhart (2001) fit, or (2.4 ± 1 .0) x 10'^ M© with the Merritt & 
Ferrarese (2001) fit, for a la allowed range of (0.9-3.4) x 10'' 
Mq. Given this result, the projected radius of the gravitational 
sphere of influence {rg = GM/ a^) of the black hole in Mrk 501 
is expected to be ~ Of 04-0"16. This raises the possibility that 
a genuine stellar-dynamical mass measurement could be done 
with the Hubble Space Telescope (HST), if M, is in the upper 
half of this uncertainty range. 

The main caveat to this result is that the upper end of the 
M, - correlation is rather uncertain at present. In the range 
M, > 10^ Mq, there are as yet only three galaxies with dynam- 
ical M, measurements derived from HST data. HST programs 
currently in progress wiU rectify this situation during the next 
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few years with many additional measurements, giving better 
constraints on the amount of intrinsic scatter in the M, - <ji, re- 
lation in this mass range. 

5. DISCUSSION 

The correlation between black hole mass and galaxy bulge 
luminosity can also be used to obtain black hole mass esti- 
mates. The Mrk 501 host galaxy has a total B magnitude of 
14.4 (Stickel et al. 1993). Correcting for Galactic extinction of 
Ab = 0.084 mag (Schlegel, Finkbeiner, & Davis 1998) and ap- 
plying a /T-correction of 0.16 mag (Pence 1976), the galaxy has 
Mb = -21 .6 mag for //q = 70 km s"^ Mpc"^ . From the fit to the 
M, -Lbuige correlation given by Kormendy & Gebhart (2001), 
the expected mass is M, = 6.1 x 10** Mq. While this value is 
outside our 1 a uncertainty range on M, from the M, — g^, re- 
lation, the scatter in the M, -Lbuige correlation is more than an 
order of magnitude in M, at fixed Lbuige (Kormendy & Geb- 
hart 2001). Also, the host galaxy luminosity is rather uncertain; 
the values compiled from the literature by Nilsson et al. (1999) 
show a range of ~ 0.5 mag among different authors. (We note 
that for the adopted luminosity, the galaxy lies very close to the 
mean Faber- Jackson relation for nearby ellipticals.) Thus, the 
level of disagreement between these two M, estimates is not a 
cause for concern; M87 is an outlier by a similar amount in the 
M, -Z-buige correlation. The M,-a^ relation is a much more 
reliable predictor of M, because its intrinsic scatter is much 
smaller, less than 40% (Gebhardt et al. 2000a) and possibly near 
zero (Ferrarese & Merritt 2000). 

Some previous estimates of the black hole mass in Mrk 501 
have been obtained by examination of the spectrum and vari- 
ability of its nonthermal emission. For example. Fan, Xie, & 
Bacon (1999) derive M, by combining the variability timescale 
(from which they estimate the physical size of the emitting re- 
gion) with the assumption that the 7-rays in blazars are pro- 
duced at ~ 200 Schwarzschild radii from the black hole. For 
Mrk 501 and seven other 7-ray blazars, they find central masses 
in the range (1 -7) x 10^ Mq. Rieger & Mannheim (2000) in- 
terpret the periodic behavior observed in X-ray and 7-ray light 
curves as evidence for a binary black hole in Mrk 501, and pro- 
pose a geometric model for a jet originating from the less mas- 
sive black hole to estimate M, < 10*^ Mq and (4-42) x 10^ 
Mq for the two components of the binary. Recently, Wang, 
Xue, & Wang (2001) have devised a method to estimate M, 
in blazars from the peak luminosity and peak frequency. Us- 
ing data for a large sample (but not including Mrk 501), they 
show that their method implies masses of ^ 10'~' ' Mq for low- 
frequency peaked blazars, and 10^"^ Mq for high-frequency 
peaked blazars. 

However, these methods have not been calibrated for galax- 
ies whose black hole masses have been measured dynamically, 
so it is difficult to assess their accuracy. Black hole masses be- 
low 10^ Mq in BL Lac objects would be difficult to reconcile 
with the M, -Lbuige correlation, since HST imaging has conclu- 
sively demonstrated that most BL Lac hosts are luminous el- 
lipticals (Urry et al. 2000). In addition, masses of < 10*^ Mq 
would pose a problem for unified models of radio-loud AGNs, 
in which BL Lacs are interpreted as radio galaxies oriented with 
the radio jet along our line of sight (e.g., Urry & Padovani 



1995). Nearby FR 1 radio galaxies, which would presumably 
appear as BL Lac objects if viewed pole-on, have black holes 
with masses in the range ~ (0.3-3) x 10^ Mq (e.g.. Harms et 
al. 1994; Ferrarese, Ford, & Jaffe 1996). Thus, for Mrk 501, the 
central mass derived from either fj* or Lbuige is consistent with 
the range of masses expected from unified schemes. 

The Eddington ratio of the active nucleus in Mrk 501 can be 
determined if its bolometric luminosity is known. Only indi- 
rect estimates of Lboi are possible, since the observed nonther- 
mal emission is highly beamed. One way to estimate Lboi is 
by comparison with FR I radio galaxies, the unbeamed coun- 
terparts of BL Lac objects. Bolometric luminosities for four 
nearby FR I radio galaxies (M84, M87, NGC 6251, and NGC 
4261) are available from Ho (1999), and nuclear emission-line 
measurements for these objects are listed by Shuder & Os- 
terbrock (1981) and Ho, Fihppenko, & Sargent (1997). For 
these four galaxies, the ratio L(Ha;-H[N II])/Lboi has a mean 
value of 4.4 x 10""' with a scatter of only ^ 25% over an or- 
der of magnitude range in Lboi. Mrk 501 has L(Hq;H-[N II]) 
= 1.1 x 10^1 ergs"' (Stickel etal. 1993). If itsL(HaH-[NlI])/Lboi 
ratio is similar to these local FR I galaxies, as might be ex- 
pected from unified models, then it has Lboi » 2.4 x lO''^ erg 
s"' and L/Ledd « 10"^. Mrk 501 is evidently an extremely sub- 
Eddington accretor; this strengthens previous indications that 
nearby TeV-emitting blazars are intrinsically weak AGNs (e.g., 
CavaUere & Malquori 1999). 

The correlations between M, and host galaxy properties 
make it possible to consider new tests of AGN unification sce- 
narios, by examining the distributions of black hole masses and 
Eddington ratios between different classes of AGNs. Such tests 
have been performed using the host galaxy luminosity to esti- 
mate M, (Treves et al. 2001), but a comparison based on stellar 
velocity dispersion would provide improved constraints on the 
central masses, and new probes of the demographics of AGN 
classes such as FR I and FR II radio galaxies and blazars. Also, 
it would be worthwhile to test whether the black hole masses 
differ systematically between high- and low-frequency peaked 
blazars, as advocated by Wang et al. (2001). It should be possi- 
ble to measure stellar velocity dispersions in several other low- 
redshift BL Lac objects. At recession velocities beyond 12,000 
km s"' , the Ca A8542 line will be redshifted into telluric H2O 
absorption bands and difficult to measure accurately, but the 
Mg lb spectral region will still be observable. We are beginning 
a spectroscopic survey of additional nearby BL Lac objects in 
order to address these questions. 
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Fig. 1 . — Optical spectrum of Mrk 501. The spectrum of the elliptical galaxy NGC 4278, observed during the same run, is shown for comparison. An arbitrary 
constant (different for the blue and red sides) has been added to the NGC 4278 spectrum. 




Fig. 2. — Best fits of broadened stellar templates to the Mrk 501 spectrum. The regions included in the calculation are shown with thick lines, and dashed lines 
denote the portions of the template spectra not used in the fit. On the red side, the strongest telluric absorption features remaining in the spectrum are labeled. 



